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It is believed that soon after the Planck era, spacetime should have a semi-classical 
nature. According to this idea, it is unavoidable to modify the General Relativity 
theory or look for alternative theories of gravitation. An interesting alternative found 
in the literature is to take into account two geometric counter-terms to regularize 
the divergences from the effective action. These counter-terms are responsible for a 
higher derivative metric theory of gravitation. In the present article, our main aim is 
to investigate if the process of isotropisation, which could well occur naturally in such 
a higher derivative metric theory, could generate gravitational waves. First of all, in 
order to proceed with such an investigation, we evaluate the time evolution of the 
resulting Weyl's scalars. Then, we follow with a perturbative approach considering 
that the background metric is given by the Minkowski spacetime. Our entire analysis 
is restricted to the particular Bianchi I case. 
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I. INTRODUCTION 

A semi-classical theory considers the back reaction of quantum fields in a classical 
geometric background. This approach started a long time ago with de Witt and since 
then, their consequences and applications are still under investigation (see, e.g., Ref. 

Diff'erently from the usual Einstein-Hilbert action, the predicted gravitational action 
allows differential equations with fourth order derivatives, which is called the full theory 
[it, see also jsj]. 

The full higher order theory was previously studied by Starobinsky ^ , and more recently, 
also by Shapiro, Pelinson and others jsf. In (sl only the homogeneous and isotropic spacetime 
is studied. 

Here, the full theory with four time derivatives is addressed, which was first investigated 



for general Bianchi I spacetimes in Tomita's article jo]. They found that the presence of 
an anisotropy contributes to the formation of a singularity. Schmidt has a recent and very 
interesting review on higher order gravity theories in connection to cosmology 

For general anisotropic Bianchi I homogeneous spacetimes, the full theory reduces to 
a system of nonlinear ordinary differential equations. The numerical solutions for this 
system with derivatives of fourth order in time was previously obtained by one of us and a 



a 



collaborator 

The question of stability of Minkowski space in this quadratic theory of gravity was 
also addressed and we have obtained that there are many initial conditions that result in 
Minkowski space asymptotically (q]. The analysis is restricted to Bianchi I type spacetimes. 
There is a basin of attraction to the Minkowski spacetime. Isotropisation can also occur 
towards de Sitter spacetime for a non zero cosmo logical constant. In this sense the above 
theory is structurally stable, since for a set of non zero measure of initial anisotropic 
geometries, isotropisation occurs. 

It should be mentioned that the isotropisation process for non zero cosmological constant 
is not a peculiarity of these higher derivatives theories. It also occurs in the Einstein's theory 



of General Relativity, as proved in a very interesting article 10|], where for A > all Bianchi 
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models, except a highly positively curved Bianchi IX, the theory becomes asymptotically a 
de Sitter spacetime. The remarkable difference between the theory of General Relativity and 
the higher derivatives theories is that the isotropisation depends on the initial conditions. 
The occurrence of isotropisation is strongly dependent on the initial conditions for quadratic 
curvature theories such as those studied here. 

In the present work our aim is to investigate how isotropisation is approached, and if 
this very process could naturally generate gravitational waves. For this reason an initial 
condition is chosen in such a way that it asymptotically generates the Minkowski spacetime, 
and the Weyl's scalars are also obtained for this solution. It must be emphasized that this 
solution, although numeric, is an exact solution of this particular quadratic gravity in the 
sense that it depends only on the precision of the machine. We then consider a perturbative 
approach to investigate if gravitational waves are in fact generated during the isotropisation 
process. 

The article is organized as follows. In section II a brief discussion on the perturbative 
approach is given. The field tensor is obtained to first order in the perturbation parameter, 
resulting in the linearized field equation. In section III we evaluate the Weyl invariants 
^'2 and \E'4 and check their time dependence, albeit having a decreasing behavior they are 
non null which indicate the presence of gravitational waves. Then, two distinct methods 
are applied to extract the GW amplitudes from the exact numerical solution. It is also 
shown that the metric perturbation can be understood in the TT gauge as a superposition 
of the Minkowski space, a pure tensorial perturbation and a pure scalar perturbation. Our 
conclusions and final remarks are presented in section IV. 

The following conventions and units are taken Rl^^ = T^^^ — Rat = Racbi ^ — 
metric signature — h ++, Latin symbols run from — 3, Greek symbols run from 1 — 3 and 



II. PERTURBATIVE ANALYSIS 



A. The full equations 



The Lagrangian density we consider here is 




(1) 
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By varying the action I = J Cd^x with respect to the metric we find the field equations 



Gn 



(1) 

ab 



aH. 



(2) 
ab 



T, 



ab 



^^9ab, 



where 



H. 



(1) 

ab 
(2) 



Gab 



^dabR — "^RRab + '^R;ab — '^'^RQabi 

■::gabRmnR'^^ + R;ab — '^R'^"' Rcbna — ORab — -ORgab, 
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Rab 



1 

^dabR, 



(2) 

(3) 
(4) 
(5) 



and Tab is the energy momentum source, which comes from the classical part of the 
Lagrangian Cc- Only vacuum solutions Tab = = will be considered in this article since 
it seems the most natural condition soon after the Planck era. The covariant divergence 
of the above tensors are identically zero due to their variational definition. The following 
Bianchi Type I line element is considered 



ds'^ = -dt^ + [ai{t)fdx'^ + [a2{t)fdy'^ + [a3{t)fdz'^ 



(6) 



which is a general spatially fiat and anisotropic space with proper time t. With this line 
element all the tensors which enter the expressions are diagonal. The substitution of Eq. 
([6]) in Eq. ([2]) with Tab = results, for the spatial part of (l2l), in differential equations of 
the type 

f_ 

dt^' 
d^ 
dt^' 
d' 



df^ 



ai = fi 



d' „ 



df^ 



as = h 



dt^ 



di^ flj, di^ di 



(7) 
(8) 
(9) 



where the functions involve the ai, 02, and their derivatives in a polynomial fashion 
(see the Appendix of the Ref. js]). The Ref. 11 1 shows that the theory which follows 



from Eq. ([T]) has a well posed initial value problem. This issue has some similarities with 
General Relativity theory since it is necessary to solve the pro blem of boundary conditions 
and initial conditions simultaneously. According to Noakes HI], once a consistent boundary 
condition is chosen initially, the time evolution of the system is uniquely specified. In Hj 
the differential equations for the metric are written in a form suitable for the application of 



the theorem of Leray 



12|. 
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For homogeneous spaces the differential equations derived from Eq. ([2]) reduce to a set 
of non hnear ordinary differential equations. Then, instead of going through the general 



construction given in Ref. [Uj, in this particular case, the existence and uniqueness of the 
solutions of Eqs. ©-(IQ]) reduce to the well known problem of existence and uniqueness of 
solutions of ordinary differential equations. For a proof on local existence and uniqueness of 
solutions of differential equations see, for instance jisl . 

Besides the Eqs. ([7])- ([9]), we have the temporal part of Eq. ([2]). To understand the role 
of this equation we have first to study the covariant divergence of the Eq. ([2]), 

Va^"'' = daE"^ + Vl^E''^ + r^^E'^" = 0. (10) 

Recall that the coordinates that we are using are those specified by the metric ([6]) and 
that Tab = 0. Since the differential equations ©-(E]) are solved numerically = 0, 

doE'^ + r:oi?°° + r°oEO° = 0, (11) 

= 0, (12) 

where £"00 is the 00 component of Eq. ([2]). If £"00 = initially, it will remain zero at any 
instant. Therefore the equation £00 acts as a constraint on the initial conditions and we 
use it to test the accuracy of our results. The full numerical solution to this problem was 
worked out in the Ref. [91]. 

B. Perturbative Approximation 

Following MTW |14] we consider metric perturbations in the form 

9ab = gif^ + ehab (13) 

g'^b ^ g<^b^-eh^^ + e^h''^hj; + ... (14) 

In a free falling frame we have g^J^^ = rjab for the background metric and locally T^^^l^ = 
so that the connection compatible with gab is identical to F^^ = eS'^ab then 



V.n = Ta.,b = Ta\b - eSl,T, (15) 
1 

2' 



Sab = -^dls) {hsa\b + hsb\a " hab\s) , (16) 
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where | is the covariant derivative with respect to the background metric g^f^ and Va =; a 
is the covariant derivative with respect to gab- The tensors which appear in Eqs. ([2]), ©-([S]) 
can be expanded in powers of e to first order, namely 

= - l^bR^^^ - i^lf i?^ (17) 
Hil^' = 2RU - 2n(^)i?^,if (18) 
hT = RU-°'''RI^-1°'''R'9[T. (19) 

where R^, i?^^, are the Ricci scalar and the Ricci tensor to first order in e, Rlfj^d is the 
Riemann tensor to first order (see [l^]) and D^T" = ^^'""V^V^T'^ = Tj'J^ 

These expressions are combined to give a perturbative sum for Eat, to first order 

Eab = E^^^ + eElf, + .. 

EL = + (/3 - Hil^' + aHj^J' - hh^b (20) 

where E^^ correspond to the background equations which is supposed to be exactly 
E^f^ = 0. By imposing E^f^ = it is possible to obtain linearized partial differential equations 
for hab- 



III. NUMERICAL RESULTS 



In what follows we apply two distinct methods in order to obtain a GW solution as a final 
result of isotropisation. Our initial ansatz is that after some arbitrary time interval we can 
consider the full solution as asymptotically Minkowskian. With this is mind, it is possible 
to split the asymptotic solution in a fiat background metric plus a perturbative quantity. 
Thus, in the first approach we numerically extract the metric perturbation from the full 
solution, while in the second approach we explicitly write down the governing equations for 
the perturbations. 



A. Decomposition of the asymptotic solution 

First consider the full numerical solution of the Eqs. (JTj), ([H]) and ([9]). In FIG. [H we 
show a particular solution for which isotropisation occurs, the parameters chosen are a = 2, 
(3 = —1, A = with initial conditions di = 0.1, 02 = 0.3, ai = 1, 02 = 1, 03 = 1; ds is 
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FIG. 1: The numerical solution of Eqs. ([7]) ([8]) ([9]) for ai{t), a2{t) and 03(4) together with a 
linear regression. The parameters are a = 2.0, (3 = —1.0, A = and the initial condition is 
hi = 0.1, (12 = 0.3, oi = 1, 02 = 1, as = 1. The numerical value of 03 is chosen in accordance with 
(jl2p . and all the other higher derivatives are chosen to be zero. The linear regression is done for 
a time interval around t ~ 8 x 10'' in natural units, and the resulting values for the coefficients 
are, see <^ bi = 0.0002293455, 62 = 0.0007793157, 63 = 0.0005822617, ci = 9173.663, C2 = 
31171.99, C3 = 23289.99. 

chosen in accordance with Eq. f ll2p . and all the other higher order derivatives are chosen to 
be zero. 

We can check the existence of gravitational waves in that solution by evaluating the 
resulting Weyl's invariant that has spin 2 (see, e.g., [l5|). The time dependency of \l/4 is 
displayed in FIG. |2l where, for completeness, we are also showing the ^2 invariant that has 
spin 0. They were evaluated for the initial condition and parameters considered above and, 
as can be seen, both of them decay with time albeit being non zero. For the particular line 
element chosen here we have \l/o = and ^1 = \E'3 = 0. 

Since the existence of tensor gravitational waves is evidenced by the presence of a non- 
null \l/4 invariant, and a scalar perturbation is also present as shown by the invariant ^21 
the next step is to find the time evolution of the amplitudes of the metric perturbations. 

In a first approach we consider an arbitrary time interval around t ~ 8 x 10'^ in which we 
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FIG. 2: Time decaying behaviour of the Weyl's invariants ^'4 and ^'2^ respectively of spin 2 and 
0, for the full solution for ai{t), 02 (t) and 03 (t) shown in FIG. [TJ 

have applied a linear regression (see FIG. [1]) and the following asymptotic expressions were 
obtained 



01 ~ Ci + felt 

02 ~ C2 + b2t 

03 ^ C3 + bst, 



(21) 



where fei = 0.0002293455, fe2 = 0.0007793157, feg = 0.0005822617, Ci = 9173.66, C2 = 
31171.99 and C3 = 23289.99. 

After the subtraction of Eq. ( 12T]) from the exact numerical solution of the Eqs. ([7]), (jH]) 
and ([9]), it can be seen that the metric of the complete solution can be rearranged in the 
following way 

9ab = ~9T + hah. (22) 

where ^^f^ is a constant background given by ^^f^ = diag(— 1, c^, Cg, Cg) and hab is a 
superposition of the time oscillating functions hi{t), h2{t) and h^{t) plus a term that gives 
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FIG. 3: The oscilating part hi{t) = [aj(t)]^ — [q + bit]'^. The values of ai are the numerical solutions 
of Eqs. ©-(111). The values of q and bi are shown in FIG. [H Note that the time origin was changed 
to t = 8.0001 X lO'^ + 

a quadratic in time contribution 

/o 

hi{t) + 2cibit + bjt'^ 
hit) + 2c2b2t + blt^ 



habit) 

















• (23) 



hit) + 2csbst + bjt^ / 



The oscillating behavior of the functions hit), hit) and hit) are shown in FIG. 
[HI Therefore, we can say that the sum + hat is identical to the metric Qab = 
diag(— 1, ai(t)^, a2(t)^, a3(t)^). Now, we must ask ourselves if the perturbation hb indeed 
represents gravitational waves. We see that hb is spatially divergenceless since hb depends 
only on time. Also, its time components are null as can be seen from the matrix fl23l) . but it 
is not traceless. Thus, in order to obtain a gravitational wave signal in TT gauge we need 
to subtract the trace of hb generating the new quantity 

1 



h 



ab 



hab - ^^if 



(24) 



which finally has the properties 



^ab 



/i = 0, 



(25) 
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FIG. 4: This is the plot of the spatial trace g^^^hij. It can be seen that the trace of h is almost 
constant in comparison with the individual amplitudes of the components of hat- Note that the 
scalar component is smaller than the tensor component. Here also a new time origin was chosen 
to be t = 8.0 X 10'^ + t+. 

and therefore we conclude that hat represents tensor gravitational waves. The trace 
h = g'^Q-^hah represents a scalar perturbation (related to ^2) that evolves smoothly with 
time and is smaller than the tensor perturbation hab as can be seen in FIG. HI 

B. The linearized solution 



Now, let us formally write the linearized system of equations which is asymptotically 

consistent with the full numerical solution. Considering the metric as gat = .^ah + hat, the 

linear system to be solved is obtained by setting E^^ = 0, resulting in 

r/4 „ f/3 . " A 

= H it) —h, + /J(t)/., + fl^{t)h, + (t)/.,- 

„ " A 

r/4 „ (j'i ^ ■■ X 

= fi, (t) —h, + fl{t)h, + fl^{t)h, + fl^{t)h, (26) 

where h is defined as the oscillatory part of the metric perturbation and the coefficients fij{t), 
also depend on the parameters a, /3, and A as can be seen in the Appendix. Now consider 
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FIG. 5: The numerical time evolution of hi for the linearized system given in Eqs. (j26p for 
an initial condition consistent with the one shown in FIG. [3l It can be seen that the solution is 
indistinguishable from FIG. El Again, note that the time origin was changed tot = 8.0001x10^+**. 

an initial condition for the linearized differential equations for hi(t). Once the numerical 
values of the and their derivatives up to the third order are initially known, the numerical 
values of the hi and their derivatives up to the third order can also be obtained. Then a 
numerical solution of the above linearized system is specified uniquely and is shown in FIG. 
[5l As can be seen, the solution of the linearized system is exactly equal to the oscillatory 
solution found in our preceding analysis, which corroborate our initial ansatz that the full 
solution can be treated as being asymptotically Minkowski. 

Again, hab cannot be regarded as tensor GWs since it contains a scalar component. Thus 
we extract its trace as we did in the preceding section and we can say that the remaining 
quantity hat = hab — 9ab^h/4 represents the GW perturbations. 
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IV. CONCLUSIONS 

In the present article it is considered general anisotropic Bianchi I homogeneous 
spacetimes given by the line element given by Eq. ([6]). For this line element, the theory given 
in Eq. ([T]), reduces to a system of ordinary nonlinear differential equations with four time 
derivatives. This system is numerically investigated. The assumption of a homogeneous 
Universe is artificial, and still presents a first generalization for its very primordial stages. 

There is a well known conjecture that dissipative processes can take place in an infinite 
dynamical system reducing the number of degrees of freedom to just a few. In some sense, 
this argument reduces the artificiality of the assumption of a homogeneous Universe. 

Only the vacuum energy momentum classical source is considered. It should be valid soon 
after the Planck era in which vacuum classical source seems the most natural condition. 

We have found that the Weyl's invariants \l/4 and \l/2 are non-null for the full numerical 
solution, although they have a decaying behavior. This is an indication of the existence 
of gravitational waves due to an initially anisotropic Universe. Then, we have adopted 
two distinct approaches in order to obtain a description for the linearized gravitational wave 
tensor hab- The two methods have coincided in the prediction of the oscillatory pattern of the 
tensor component of gravitational waves. A scalar component is also present, that evolves 
smoothly in the time interval considered and is smaller than the tensor component. We 
conclude that the asymptotic geometry of an initially anisotropic spacetime is a superposition 
of Minkowski space and gravitational waves with tensor and scalar components. 
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Appendix A: Linearized System 



Recall that the perturbation is hnearly related to the metric while the scale factors are 
quadratic, then 

hi^aj- (ci + bitf 
h^^al- (c2 + b2tf 

^3 = 03 - (C3 + 63^)^ 

hi — 2aiai — 26i(ci + bit) 
h2 = 20202 - 262 (C2 + b2t) 
ha = 20303 - 263(03 + b^t) 
hi = 2oiai + 2(oJ^ - 2bl 
h2 = 20202 + 2(02)^ - 262 
h3 = 20303 + 2(03)^ - 263 
-hi^2ai-ai + 6aiai 

— /i2 = 202^02 + 60202 
dt"^ dt"^ 

d^ d^ 

^/.3 = 203^03 +60303^ 

|j/., = 2o,|jO, + 8d,|^o, + 6(d,f 

(^4 ^4 ^3 

^/l2 = 202^02 + 802^02 + 6(02)^ 

|j/.3 = 203^,03 + 8d3|,03 + 6(d3r. 



Now the differential equations given in Eq. (126|1 



yi 




y2 


=k 


y3 








2/4 


= d^^' 


2/5 


=k 


2/6 


=k 


2/7 


=k 






2/8 




2/9 


= ^3 


2/10 


= ^3 


2/11 


= ^3 
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—hi = (12 /3 A ci^C2^yg - a A ci^C2^|/9 - a A C2^cz'yi + 12 A ci^ca^ys/? - A cici^y^a- 

24 ^ A C2^C3^yi - 2 q; A ca^cg^ftiCii - 2 a A ci^ca^fcsCs^ + 24 /3 A ci^C2%zCzt- 
48 ^ A ca^cg^fcicit + 24 A ci^C3262C2^/5 - 2 A Cri^h2C2ta + Ci2c2^?/iiQ;+ 
6 /3 ci'cs'i/T + a ci'ca'yr - 12 ca'cs'i/s/? + 12 ci^C2%^^P + 2 cs'cg'fei'a- 
24 C2^cz%i^^ + cs^cs^yao; + 2 a 01^03^62^ + 2 01^02^63^0; + 12/3 01^03^62 V 
6ci^C2Vi/5 - q;Aci^C2^63V - q;Ac2^C3^6iV + 12Aci^C3%2h'^p- 
24/3Ac2^C3^6iV + 12/3Aci2c2^63V - Aci2c3^62Va) / (l8G'c3Va/3) 

^/t2 = (12 ^ A ci^C2^?/9 - a A ci^C2^y9 - a A C2^c^^yi - 24 A ci^cs^ysP - A 01^03^^50;+ 
12 /3 A C2^C3^|/i - 2 a A cs^Cg^feiCit - 2 a A Ci^Cs^^sCgt + 24 /3 A 01^02^6303^+ 
24^Ac2^C3^6iCit - 48 A 01^03^6202^/5 - 2 Aoi2o326202to + 01^02^10;- 

12 p 0i'03'?/7 + a 0i'03'?/7 + 6 02'03'?/3/5 + 12 Ci''c2%^P + 2 02'03'6i'q;+ 

12 02^03^61^/3 + 02^03^i/3« + 2 a 01^03^62^ + 2 01^02^63^0- 
24 /3 ciWb2^ + 6 ci^C2^yul3 - a A oi2o2'63'i' - a A 02'o3'6i2i2_ 
24A01V62V/3 + 12/3Ao2^03%V + 12^AoiW63V- 
Aci'o3'62't'a) / (l8Goi2o3'a;g) 



-^ha = (-24 P A oi^os^yg -aA ciC2^yg -ah 02^03^yi + 12 A oi^og^ys/^ - A 01^03^^50;+ 
12 /3 A 02^03^1/1 - 2 a A 03^03^6101^ - 2 a A Oi^02^6303t - 48 /3 A 01^02^6303^+ 
24 /3 A 02^03^6iOii + 24 A 01^03^6202^/? - 2 A 01^03^6202^ + 01^02^1/110+ 
6 ,5 ci'c3'y7 + a Ci'c3'y7 + 6 C2^C3^y3p - 24 Ci'c2'63'^ + 2 C2'c3'6i'o;+ 
12 02^03^61^/3 + 02^03^i/3a + 2 a 01^03^63^ + 2 01^02^63^0; + 12/3 01^03^62=^- 
12 01^02^1/3 - a A 01^02^63 V - a A 02^03^61 V + 12 A 01^03262^^^/3+ 
12^ Ac2'c3'6iV - 24^ Aci2c2'63'i' - A 01203^62^0;) / (18^0220120;^) 
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